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Properties of Cold Dark Matter Halos at z > 6 

Renyue Cen^, Feng Dong^, Paul Bode^, and Jeremiah P. Ostriker^ 

ABSTRACT 

We compute the properties of dark matter halos with mass 10^'^ — 10^ Mq 
at redshift 2; = 6 — 11 in the standard cold dark matter cosmological model, 
utilizing a very high resolution N-body simulation. We find that dark matter 
halos in these mass and redshift ranges are significantly biased over matter with 
a bias factor in the range 2 — 6. The dark matter halo mass function displays a 
slope of 2.05 ±0.15 at the small mass end. We do not find a universal dark matter 
density profile. Instead, we find a significant dependence of the central density 
profile of dark matter halos on halo mass and epoch with = 0.4 — 1.0; the 
high-mass (M > 10^ M0) low-redshift {z ~ 6) halos occupy the high end of the 
range and low-mass (M ~ 10'^ M©) high-redshift {z ~ 11) halos occupy the low 
end. Additionally, for fixed mass and epoch there is a significant dispersion in ao 
due to the stochastic assembly of halos. Our results fit a relationship of the form 
ao = 0.75((1 + 2)/7.0)-^-2^(M/10^Mq)°-i^(i+^)/^-° with a dispersion about this fit 
of ±0.5 and no systematic dependence of variance correlated with environment. 
The median spin parameter of dark matter halos is 0.03 — 0.04 but with a large 
lognormal dispersion of ~ 0.4. Various quantities are tabulated or fitted with 
empirical formulae. 

Subject headings: cosmology: theory — intergalactic medium — large-scale struc- 
ture of universe — quasars: absorption lines 
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1. Introduction 

The reionization epoch is now within the direct observational reach thanks to rapid 
recent observational advances in two fronts — optical quasar absorption from Sloan Digi- 
tal Sky Survey (SDSS) (Fan et al. 2001; Becker et al. 2001) and the Wilkinson Microwave 
Anisotropy Probe (WMAP) experiment (Kogut et al. 2003). The picture painted by the 
combined observations, perhaps not too surprisingly, strongly suggests a complex cosmologi- 
cal reionization process, consistent with the double reionization scenario (Can 2003). It may 
be that this is the beginning of a paradigm shift in our focus on the high redshift universe: 
the star formation history of the early universe can now be observationally constrained. 

It thus becomes urgent to theoretically explore galaxy and star formation process at 
high redshift in the dark age {z > 6). In the context of the standard cold dark matter 
model it is expected that stars within halos of mass 10^ — 10^ at high redshift play an 
important, if not dominant, role in determining how and when the universe was reionized. 
Furthermore, these fossil halos may be seen in the local universe as satellites of giant galaxies. 
This linkage may potentially provide a great leverage to nail down the properties of the high 
redshift galaxies. 

In this paper, as a step towards understanding galaxy formation at high redshift, we 
investigate the properties of dark matter halos at 2; > 6, using very high resolution TPM 
N-body (Bode et al. 2001; Bode & Ostriker 2003) simulations. While there is an extensive 
literature on properties of halos at low redshift, there is virtually no systematic study of 
dark halos ai z > 6. The LCDM simulation has a comoving box size of 4/i~^Mpc with 
512^ = 10^-^ particles, a particle mass of nip — 3.6 x 10^ Mq, and comoving gravitational 
softening length of 0.14 h~^kpc. These resolutions allow us to accurately characterize the 
properties of halos down to a mass 10^'^ Mq (having about 100 particles within the virial 
radius). The outline of this paper is as follows. The simulation details are given in §2. In §3 
we quantify properties of dark matter halos in the mass range 10^'^ — 10^ Mq, including the 
mass function, bias and clustering properties, density profile distribution, angular momentum 
spin parameter distribution, internal angular momentum distribution and peculiar velocity 
distribution. We conclude in §4. 



2. The Simulation 

A standard spatially fiat LCDM cosmology was chosen, with Q^n = 0.27 and Q\ = 0.73; 
the Hubble constant was taken to be 70 km/s/Mpc. The initial conditions were created using 
the GRAFIC2 package by Bertschinger (2001). The matter transfer function was calculated 
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with the included Boltzmann integrator (Ma & Bertschinger 1995), using Qfe/i^ = 0.211 for 
the baryon fraction and ag — 0.73 for the normahzation of the matter power spectrum. 

The simulation contained N — 512^ particles in a comoving periodic box Ah'^Mpc 
on a side, making the particle mass nip — 3.57 X 10^ H-^Mq. The starting redshift was 
2; = 53, and the system was evolved down to 2; = 6. The evolution was carried out with 
the parallel Tre(^Particlc-Mcsh code TPM (Xu 1995; Bode, Ostriker, & Xu 2000; Bode & 
Ostriker 2003), using a 1024"^ mesh. The evolution took 1150 PM steps, with particles in 
dense regions taking up to 19,500 steps. The run was carried out using up to 256 processors 
on the Terascale Computing System at Pittsburgh Supercomputing Center. 

The spline softening length in the tree portion of the code was set to e = 0.14 /i~^kpc. 
With this softening length, relaxation hj z = 6 inside the core of a collapsed halo (assuming 
an NEW density distribution with c = 12) will not be significant over the course of the 
simulation for those objects containing more than 100 particles. The opening angle in the 
Barnes-Hut criterion used by TPM was 9 — 0.577, and the time step parameter rj — 0.3; also, 
the initial values for locating trees were A — 2.0 and B — 12.5 — see Bode & Ostriker (2003) 
for details. In the TPM code, not all regions are treated at full resolution. The limiting 
density (above which all cells are put into trees for increased resolution) rises with time. By 
the end of this run, all cells containing more than 18 particles are still being followed at full 
resolution. Thus this factor is not important if the analysis is limited to halos with over 100 
particles. 

Dark matter halos are identified using DENMAX scheme (Bertschinger & Gelb 1991), 
smoothing the density field with a Gaussian length of 300/i~^kpc. In computing all quantities 
we include all particles located inside the virial radius of a halo. 

3. Results 
3.1. Pictures 

Eirst, we present visually a distribution of the dark matter mass and dark matter halos 
of varying masses. Eigure 1 shows the distribution of dark matter particles projected onto 
the x-y plane. Eigures (2a,b,c,d,e) show the distributions of dark matter halos with masses 
greater than (10®, 10®'^, 10^, 10^'^, 10*^) M©, respectively, at 2; = 6. The progressively stronger 
clustering of more massive halos is clearly visible in the display but we will return to the 
clustering properties more quantitatively in §3.3. It is also noted that voids are progressively 
more visible in the higher mass halos than in low mass halos. 
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Fig. 1. — The distribution of dark matter particles at 2; = 6 projected onto the x-y plane 
(0.25% of the total). 
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2. — The distributions of all dark matter halos with masses greater than 
106•^ 10^ 10^•^ 10^) M0, respectively, at ^ = 6. 



-6- 



M > lO^-^ h-^ M, 




12 3 4 



Fig. 2. — Continued. 



- 7- 



M > 10^ h-'^MQ 



/tV-^;. V . ■ ■, V '.-. - - 




Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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3.2. Dark Matter Halo Mass Function 

Figure 3 shows the halo mass functions at four redshifts. Table 1 summarizes the fitting 
parameters for a Schechter function of the form 

n{M)dM = rzo(^)-" exp (-M/M.)^. (1) 

We see that the Schechter function provides a good fit to the computed halo mass function. 
The faint end slope is a = 2.05 ±0.15, consistent with the expectation from Press-Schechter 
theory (Press & Schechter 1974). While there appears to be a slight steepening of the slope 
at the low mass end from —1.9 to —2.2 from z = 6 to 2; = 11.1, it is unclear, however, 
how significant this trend is, given the adopted, somewhat degenerate fitting formula. The 
turnover at Mh ~ 10^-^ indicates the loss of validity of our simulation at the low mass 
end. 



3.3. Bias of Dark Matter Halos 

We characterize the relative distribution of halos over the total dark matter distribution 
by the following relation: 

-J^I^ = 6(M,z)(— ^)^(^'^), (2) 

<nh> < Pm> 

where Uh and < Uh > are the halo density and mean halo density; pm and < pm > are 
the mass density and mean mass density; c{M,z) is fixed to be unity at ^^'"^ > 1. This 
empirical fitting formula is motivated by the found result that there appears to be a break in 
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^^^^ at ^p""^ ~ 1- Tables 2 and 3 list the parameters b{M,z) and c{M,z). The smoothing 
length used here is 0.3 Mpc. Figure 4 shows four typical cases to indicate the goodness 
of the fitting formula. At < 1 our results (Table 2) indicate oc (^^)^"^, a 

° <pm> ^ ' <nh> ^<pm>' ' 

Table 1. Halo Mass Function Fitting Parameters 



Parameters 


z=6.0 


z=7.4 


z=9.08 


z=11.096 


no {h? Mpc-3) 


0.85 


1.20 


1.75 


2.40 


a 


1.9 


2.0 


2.1 


2.2 


{h-^ Mpc) 


8 X 10* 


6 X 10*^ 


4 X 10* 


2 X 10* 
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Fig. 3. — The halo mass functions at four redshifts. Dotted hues represent the fitted 
Schechter functions with parameters summarized in Table 1. 
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Fig. 4. — Bias of dark matter halos in four randomly selected cases upon different mass and 
redshift. The darkened lines show the fitting formulae of the bias relation. 
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Table 2. Halo Bias : b(M,z) 



Halo Mass {h'^ Mq) z=6.0 z=7.4 z=9.08 z=11.096 



>106° 1.16±0.12 1.36±0.18 1.70±0.20 2.46±0.23 

>10^-^ 1.21±0.13 1.45±0.17 2.09±0.30 3.05±0.33 

>10^° 1.37±0.19 1.60±0.23 2.77±0.40 3.92±0.90 

>10^-^ 1.72±0.32 2.12±0.47 3.46±0.87 6.17±3.17 

>10S° 2.08±0.55 2.58±0.90 4.43±1.47 9.90±4.58 



Table 3. Halo Bias : c(M,z) 



Halo Mass {h'^ Mq) z=6.0 z=7.4 z=9.08 z=11.096 



>106-° 3.00±0.55 3.80±0.79 4.14±0.80 4.90±1.49 

>10^-^ 3.27±0.71 4.20±1.04 4.41±1.31 5.52±2.43 

>10^° 3.43±0.70 4.37±1.07 4.97±2.59 7.33±6.29 

>10^-^ 3.93±1.70 4.98±2.96 5.09±3.31 7.41±9.09 

>10^° 4.18±2.28 5.61±4.68 6.32±6.12 6.77±8.90 
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b'=lhh/fm„. as fn of mass and z 
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Fig. 5. — shows the ratio of dark matter halo correlation function over the dark matter mass 
correlation function at four redshifts, z = 6.0 (top right panel), z = 7.4 (top left), z = 9.1 
(bottom right) and 2; = 11 (bottom left). In each panel six curves are shown for halos more 
massive than 10^ M© 10^-^ Mq 10^ M© W'^Mq 10^ Mq and 10^'^ Mq, respectively. 
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Fig. 6. — shows the bias as a function of halo mass at redshifts z = (6,7.4,9.1,11.1) at 
separation of lh~^Mpc as symbols (x). The curves are computed using the analytic method 
of Mo & White (1996). 



-16- 



rather rapid drop. As expected, the drop-off is more dramatic for larger halos, as visible in 
Figure 2. This implies that at ^ > 6 halos are unlikely to be found in underdense regions (on 
a scale of ~ 0.3Mpc/h). The increase of c{M, z) with redshift implies that voids are emptier 
at higher redshifts. 

Another way to characterize the relative distribution of dark matter halos over mass is to 
compute the ratio of the correlation functions, which are shown in Figure 5. The correlation 
function ^(r) is calculated by counting the number of pairs of either particles or halos at 
separation r (using logarithmically spaced bins) and comparing that number to a Poisson 
distribution. It can be seen that the bias falls in the range 2 — 6, with the trend that the 
more massive halos are more biased and at a fixed mass halos at higher redshifts are more 
biased, as expected. Figure 6 recollects the information in Figure 5 and shows the bias as 
a function of halo mass at four different redshifts at the scale chosen to be l/i^^Mpc. The 
agreement between our computed results and that using analytic method of Mo & White 
(1996) is good, indicating that the latter is valid for objects at scales and redshifts of concern 
here. 



3.4. Dark Matter Halo Density Profile 

We use a variant fitting formula based on the NFW (Navarro, Frenk, & White 1997) 
profile: 



P'^ (-l—\a{lj^_r_\A-'ia' ("^^ 

Note that for NFW profile, a = 1. An important difference, however, is the scaling radius 
used. We use the radius where the logarithmic slope of the density profile is —2, r_2, instead 
of the more conventional "core" radius. This is a two parameter fitting formula, a and r_2, 
while Ps is a function of a and r_2 at a fixed redshift, since the ovcrdcnsity interior to the 
virial radius r„ is assumed to be known. This fitting formula is intended for the range in 
radius r < only. 

We fit the density profile of each halo using the least squares method. Four randomly 
selected examples of such profiles along with the fitted curve using Equation (3) arc shown in 
Figure 7, indicating reasonable fits in all cases. Both fitting parameters, a and r_2, however, 
display broad distributions. We found that there is only weak correlation between a and 
r_2. Figure 8 shows histograms for the distributions of a, for four typical cases. We fit a 
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Fig. 7. — Density profiles of four randomly selected halos over a mass range at z = 6. The 
dotted lines represent the best-fit modified NFW relation given in Section 3. 



z = 6.0 
M^: 10'° - 10'= h-' Mj, 

a„ = 0.808 ± 0.008 
a„ = 0.581 ± 0.006 




z = 6.0 
M^: 10" - 10»" h-i 

a„ = 0.896 ± 0.011 
CT. = 0.432 ± 0.008 




z = 11.096 
M^: 10'° - 10'= h-' Mj, 

a„ = 0.415 ± 0.014 
CT = 0.510 ± 0.010 




z = 11.096 
M^: 10" - 10»° h-i 

= 0.497 ± 0.037 
a = 0.558 ± 0.034 




Fig. 8. — The distributions of inner slope parameter a in the halo density profile fitting for 
four randomly selected cases upon different mass and redshift. The Gaussian fits are shown 
as smooth curves. 
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distributions using a Gaussian distribution function: 

The Gaussian fits arc shown as smooth curves in Figure 8, demonstrating that the proposed 
Gaussian fits arc good. Tables 4,5 list fitting parameters aQ and a^o, respectively. We 
recollect the data in Table 4 and shows in Figure 9 the median inner density slope as a 
function of dark matter halo different mass at four different redshifts (symbols) . The curves 
in Figure 9 are empirical fits using the following formula 

ao = 0.75((1 + ^)/7.O)-^-25(M/lO^M0)°-ii(i+^)/^-°. (5) 

It is seen that this fitting formula provides a reasonable fit for the simulated halos. 

Figure 10 shows histograms for the distributions of r_2 for four typical cases. We fit r_2 
distributions using a lognormal function: 

^(^■-2) = 7= exp (-^ (6) 

which are seen to provide reasonable fits to the data in Figure 10. Tables 6,7 fist fitting 
parameters r°2 and cr^.j, respectively. We find no visible correlation between a and r_2, as 
shown in Figure 11. 

An important point to note is that the average slope of the density profiles of small 
halos ranges from 1.0 to 0.4 from 2; = 6 to 2; = 11 (Table 4 and Figure 8), which is somewhat 
shallower than the universal density profile found by Navarro et al. (1997). Our results 
at 2; = 9 — 11 are in good agreement with Ricotti (2003), who simulated a l/i~^Mpc box 
for small halos ai z = 10. If the theoretical argument for the dependence of halo density 
profile on the slope of the initial density fiuctuation power spectrum (Syer & White 1988; 
Subramanian, Gen, & Ostriker 2000) is correct, as adopted by Ricotti (2003) to explain the 



Table 4. Density Profile : ao{M, z) 



Halo Mass {hr^ M©) z=6.0 z=7.4 z=9.08 z=11.096 

10^ ° - 10^-5 O.SliO.Ol 0.63±0.01 0.50±0.01 0.42±0.01 

2Q7.5 _ ^gs o 0.90±0.01 0.75±0.02 0.62±0.02 0.50±0.04 

> lO^ o 1.09±0.02 0.91±0.02 0.77±0.03 0.65±0.06 
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Fig. 9. — shows the median inner density slope as a function of dark matter halo different 
mass at four different redshifts (symbols). The curves are fits (Equation 5). 



Table 5. Density Profile : (Ja(M, z) 



Halo Mass {h'^ Mq) 


z=6.0 


z=7.4 


z=9.08 


z=11.096 


107.0 . 


0.58±0.01 


0.54±0.01 


0.52±0.01 


0.51±0.01 


107.5 . 108.0 


0.43±0.01 


0.45±0.01 


0.46±0.02 


0.56±0.03 




O.lliO.Ol 


0.11±0.02 


0.11±0.0.3 


0..38±0.06 
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Fig. 10. — The distributions of radius parameter r_2 in the halo density profile fitting for 
four randomly selected cases upon different mass and redshift. The lognormal fits are shown 
as smooth curves. 
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dependence of inner density slope on halo mass, it then follows that the neglect of density 
fluctuations on scales larger than his box size of l/t~^Mpc in Ricotti (2003) would make the 
density profiles of the halos in his simulation somewhat shallower than they should be in the 
CDM model. Thus, the difference in the box size (1/i^^Mpc in Ricotti 2003 versus 4/i^^Mpc 
for our simulation box) would have expected to result in a slightly steeper inner slope in our 
simulation, which is indeed the case. 

Another point to note, which is not new but not widely known, is that there is a 
large dispersion in the inner slope of order 0.5 due to the intrinsically stochastic nature of 
halo assembly. This was found earlier by Subramanian et al. (2000). Therefore, while a 
"universal" profile is informative in characterizing the mode, a dispersion would be needed 
to give a full account. This is particularly important for applications where the dependence 
on the inner slope is very strong, e.g., strong gravitational lensing. More relevant for our case 
of small halos is that, for example, the fraction of small halos with inner slope close to zero 
(i.e., flat core) is non- negligible at 2; = 6. A proper statistical comparison with observations 
of local dwarf galaxies, however, is not possible with the current simulation without evolving 
small galaxies to 2; = 0. 

Also worthwhile is to understand whether or not there is some dependence of the inner 
slope of the density profile on the central density of a halo or the environmental density where 
a halo sits. In Figures 12 and 13 we show the correlation between the inner slope of the 
density profile, a, and the central density of the halo, and between a and the environmental 
density, respectively. We find no visible correlations between either pair of quantities. But a 
relationship between halo shape and environment might have been missed due to stochastic 
variations. Thus, we have checked to see if deviation Aa, between computed and predicted 
(based on mass and epoch using equation 5) slope exists. Figure 14 shows this and no 
correlation is seen. We note that the abundances of halos in the mass range of interest 
here are on the rise in the redshift range considered (Lacey & Cole 1993), as evident in 
Figure 3. During this period halos considered the inner density profiles of halos steepen with 



Table 6. Density Profile : r°2(M,z) 



Halo Mass {hr^ Mq) 


z=6.0 


z=7.4 


z=9.08 


z=11.096 


107.0 . 107-5 


0.35±0.01 


0.40±0.01 


0.43±0.01 


0.44±0.01 


107.5 . 108.0 


0.35±0.01 


0.38±0.01 


0.39±0.01 


0.37±0.01 


> lO^o 


0.34±0.01 


0.35±0.01 


0.34±0.01 


0.34±0.02 
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Fig. 11. — A scatter plot of the inner slope of the density profile, a, versus r_2 for all the 
halos with M > 10^'^ H'^Mq at z = 6. 
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Fig. 12. — A scatter plot of the inner slope of the density profile, a, versus the halo central 
density, defined as the density at r < 0.2r^. for all the halos with M > 10^ h~^MQ at ^ = 6. 
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Fig. 13. — A scatter plot of the inner slope of the density profile, a, versus the environmen- 
tal density, defined as the dark matter density smoothed by a gaussian window of radius 
0.3/i-^Mpc, for all the halos with M > 10^ H'^Mq at ^ = 6. 
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time, consistent with the increase of logarithmic slope of the power spectrum with time, 
that corresponds to the evolving nonlinear mass scale. However, at some lower redshift not 
probed here, low mass halos will cease to form. Subsequently, the evolution of halo density 
profile may show distinct features and some conceivable correlations, not seen in Figures 12 
and 13, may show up. We will study this issue separately. 

Since the central density of a halo may be considered a good proxy for the formation 
redshift of the central region, the non-correlation between a and the central density indicates 
that, for halos of question here, the subsequent process of accretion of mass onto halos is 
largely a random process, independent of the density of the initial central "seed". The 
fact that halos of different masses show a comparable range of central density (not shown 
in figures) suggest that halos of varying masses form nearly simultaneously, dictated by the 
nature of the cold dark matter power spectrum at the high-k end; i.e., density fluctuations on 
those scales involved here depend weakly (logarithmically) on the mass. The non-correlation 
between the inner slope of the halo density profile and the environmental density may be 
interpreted in the following way. One may regard regions of different overdensities as local 
mini-universes of varying density parameters. The independence of the inner slope on local 
density is thus consistent with published results that the halo density profiles in universes of 
different Q^m do not significantly vary. 



3.5. Spin Pcirameter For Dark Matter Halos 

We compute the spin parameter defined as 

(Peebles 1969), where G is the gravitational constant; M is the total mass of the dark matter 
halo; J is the total angular momentum of the dark matter halo; E is the total energy of 



Table 7. Density Profile : (7r_2(M,z) 



Halo Mass {hr^ Mq) z=6.0 z=7.4 z=9.08 z=11.096 

10^-0 - 10^-5 0.34±0.01 0.32±0.01 0.29±0.01 0.27±0.01 

10^-5 - 10^-° 0.31±0.01 0.31±0.01 0.27±0.01 0.27±0.01 

> lO^o 0.33±0.01 0.25d=0.01 0.26±0.02 0.27±0.05 
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Fig. 14. — A scatter plot of the difference between measured inner slope of the density 
profile, a, and the fitted inner slope using equation (5), versus the environmental density. 
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z = 6.0 
M^: 10""° - 10' == h-' 

= 0433 ± 0.0004 
(Tj = 0.425 ± 0.005 
e„ = 0.0125 




M^: 10'» - 10»°» h-' Mj, 

'^medi.n = 0.0405 ± 0.0005 
<Tj = 0.422 ± 0.008 
0.0125 




z = 11.096 
M^: 10'" - 10'-25 h-' Mg 

■''m.di.n = 0.0402 ± 0.0007 
= 0.410 ± 0.011 
e. = 0.0125 




z = 11.096 
Mj: 10'«' - lOS" h-' Mg 

^„.di.„ = 0.0347 ± 0.0019 
<7j = 0.322 ± 0.034 
0.0125 




Fig. 15. — The distributions of halo spin parameter A for four randomly selected cases upon 
different mass and redshift. The modified lognormal fits are shown as smooth curves. 
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Fig. 16. — A scatter plot of the spin parameter, A, versus the halo central density, defined 
as the density at r < 0.2r„. for all the halos with M > 10^ H'^Mq at z = 6. 
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Fig. 17. — A scatter plot of the spin parameter, A, versus the environmental density, defined 
as the dark matter density smoothed by a gaussian window of radius 0.3/i~^Mpc, for all the 
halos with M > 10^ h'^Mp, at ^ = 6. 
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Fig. 18. — Comparison of the halo density profile slope parameter a versus spin parameter 
A. No correlation between the two is observed. 



-32- 



Table 8. Spin Parameter : XmedianO^,"^) 



Halo Mass (h-^ Mq) z=6.0 



z=7.4 



z=9.08 



z=11.096 



107.00 _ 107.25 
107.25 . 107.50 
107.50 _ 108.00 



0.043±0.001 
0.042±0.001 
0.041±0.001 
0.035±0.001 



0.043±0.001 
0.041±0.001 
0.039±0.001 
0.033±0.001 



0.042±0.001 
0.037±0.001 
0.035±0.001 
0.031±0.002 



0.040±0.001 
0.037±0.001 
0.035±0.002 
0.031±0.007 



Table 9. Spin Parameter : ax{M,z) 



Halo Mass (h-^ M©) z=6.0 z=7.4 z=9.08 z=11.096 



107.00 _ 


107.25 


0.43±0.01 


0.44±0.01 


0.44±0.01 


0.41±0.01 


107-25 . 


107.50 


0.43±0.01 


0.42±0.01 


0.38±0.01 


0.41±0.02 


107.50 _ 


108.00 


0.42±0.01 


0.41±0.01 


0.41±0.01 


0.32±0.03 


> 10* 


^.00 


0.40±0.01 


0.39±0.02 


0.31±0.04 


0.43±0.15 



Table 10. Angular Momentum Profile : Ao(M,z) 



Halo Mass (h'^ M©) z=6.0 z=7.4 z=9.08 z=11.096 

107.0 _ 107.5 0.51±0.03 0.49±0.04 0.44±0.06 0.44±0.08 

107.5 . 108.0 0.55±0.06 0.52±0.12 0.48±0.14 0.50±0.16 

> 10S-° 0.56±0.09 0.51±0.21 0.49±0.24 0.45±0.32 



Table 11. Angular Momentum Profile : (7^(M,z) 



Halo Mass {h~^ Mq) 


z=6.0 


z=.7.4 


z=9.08 


z=11.096 


107.0 . 107.5 


0.27±0.03 


0.28±0.04 


0.27±0.06 


0.25±0.07 


107.5 . 108.0 


0.26±0.05 


0.29d=0.10 


0.29±0.11 


0.26±0.12 


> lO^o 


0.27±0.08 


0.35±0.18 


0.29±0.16 


0.36±0.23 
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the dark matter halo; all quantities are computed within the virial radius. We fit the A 
distributions using a modified lognormal function: 

P(A) = — 4^exp(-M^±|I^i^), ^s) 

(A + eo)v27r(7A ^crx 

where Sq is fixed to be 0.0125, which is determined through experimentation. The modified 
lognormal fits for A are shown as smooth curves in Figure 15; the goodness of the fits is 
typical. Tables 8,9 lists fitting parameters Aq and (Ja, respectively. Note that the median 
value of A is \med = Aq - eo. 

We see that the typical spin parameter has a value 0.03 — 0.04. However, the distribution 
of the spin parameter among halos is very broad, with a lognormal dispersion of ~ 0.4. This 
implies that consequences that depend on the spin of a halo are likely to be widely distributed 
even at a fixed dark matter halo mass. Such consequences may include the size of a galactic 
disk and correlations between dark matter halo spin (conceivable misalignment between spin 
of gas and spin of dark matter would complicate the situation) and other quantities. 

In Figures 16 and 17 we show the correlation between A and the central density of 
the halo, and between A and the environmental density, respectively. We find that there 
may possibly exist a weak correlations between A and the central density of the halo, in 
the sense that halos with higher central densities (or equivalcntly earlier formation times) 
have lower A, with a very large scatter, whereas no correlation is discernible between A and 
the environmental density. Finally, Figure 18 show the relation between A and a, where no 
correlation is visible. 



3.6. Angular Momentum Profile For Dark Matter Halos 

Next, we compute the angular momentum profiles for individual dark matter halos. We 
then fit the each angular momentum profile by the following function in the small \j\ regime: 

M(< j) _ J M{< 0) 

where A and M(< 0) are two fitting parameters; My is the virial mass and jo = J/My = 
XGM^/y\E\^/^. 

In order to compute the angular momentum profile an appropriate smoothing window 
needs to be apphed to dark matter particles. We find that M(< 0)/Mv varies for different 
smoothing scales, with typical values around 0.2, while A remains roughly constant for each 
individual halo. However, A is broadly distributed for all dark matter halos. Figure 19 




Fig. 19. — The distributions of slope parameter A in the halo angular momentum profile 
fitting for four randomly selected cases upon different mass and redshift. The lognormal fits 
are shown as smooth curves. 
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Fig. 20. — The distributions of dark matter halo pecuhar velocity for four mass bins at 
z — 6. The lognormal fits are shown as smooth curves. 
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show histograms for the distributions of Aq, for four randomly selected cases. We fit the A 
distribution using a modified lognormal function: 



[ln{A + e^) - In A, 



^^^^ (A + eA)V2^aA^'''^'' 2(7^ ^' 

where = 0.4 is fixed through experimentation. Tables 10,11 list fitting parameters Aq and 
a A, respectively. 

In general we find our fitting formula (Equation 10) provides a good fit for each individ- 
ual halo. The distribution of matter at small j is most relevant for the formation of central 
objects, such as black holes (e.g., Colgate et al. 2003) or bulges (e.g., D'Onghia & Burkert 
2004). Our calculation indicates that the fraction of mass in a halo having specific angular 
momentum less than a certain value is roughly 0.5 times the ratio of that value over the 
average specific angular momentum of the halo. 
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3.7. Bulk Velocity of Dark Matter Halos 

Finally, we compute the peculiar velocity of dark matter halos. We find that the distribu- 
tion once again can be fitted by lognormal distributions as equation (9). In order to provide 
a good fit for the results ai z — 6 shown in Figure 20, it is found that e = 40km/s in equation 
(8) with median velocity ~ 38 ± 2km/s and lognormal dispersion (7„ = 0.22 ± 0.01. We 
caution, however, the absolute value of the peculiar velocity of each halo, unlike the quan- 
tities examined in previous subsections, may be significantly affected by large waves not 
present in our simulation box. Adding missing large waves should increase the zero point e 
to a larger value. Therefore, the peculiar velocity shown should be treated as a lower limit. 
In other words, expected peculiar velocity of dark matter halos at these redshifts are likely 
in excess of 30 — 40km/s. 



4. Conclusions 

Using a high resolution TPM N-body simulation of the standard cold dark matter 
cosmological model with a particle mass of rrip = 3.57 x 10'^ /i~^Mq and a softening length 
of e = 0.14 /i^^kpc in a 4/i^^Mpc box, we compute various properties of dark matter halos 
with mass 10®'^ — 10^ at redshift z = 6 — 11. We find the following results. 

(1) Dark matter halos at such small mass at high redshifts are already significantly biased 
over matter with a bias factor in the range 2 — 6. 
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(2) The dark matter halo mass function displays a slope at the small end 2.05 ± 0.15. 

(3) The central density profile of dark matter halos are found to be in the range (0.4 — 1.0) 
well fitted by ao = 0.75((1 + z) /7.0)-i-25(M/10^ yi^^o.m+z)/7.o ^ dispersion of ±0.5, in 
rough agreement with the theoretical arguments given in Ricotti (2003) and Subramanian 
et al. (2000). 

(4) The median spin parameter of the dark matter halos is 0.03 — 0.04 but with a lognormal 
dispersion of ~ 0.4. The angular momentum profile at the small end is approximately linear 
with the fraction of mass in a halo having specific angular momentum less than a certain 
value is roughly 0.5 times the ratio of that value over the average specific angular momentum 
of the halo. 

(5) The dark matter halos move at a typical velocity in excess of 30 — 40km/s. 

This research was supported in part by AST-0206299 and NAG5-13381. The computa- 
tions were performed on the National Science Foundation Terascale Computing System at 
the Pittsburgh Supercomputing Center. 

REFERENCES 
Becker, R.H., et al. 2001, AJ, 122, 2850 

Bcrtschinger, E.,& Gelb, J.M. 1991, Computers in Physics, 5, 164 

Bertschinger, E. 2001, ApJS, 137, 1 

Bode, P., & Ostriker, J.P. 2003, ApJS, 145, 1 

Bode, P., Ostriker, J.P., & Xu, G. 2000, ApJS, 128, 561 

Cen, R. 2003, ApJ, 591, 12 

Colgate, S.A., Cen, R., Li, H., Currier, N., & Warren, M.S. 2003, ApJL in press, astro- 
ph/0310776 

D'Onghia, E., & Burkert, A. 2004, astro-ph/0402504 

Fan, X., et al. 2001, AJ, 122, 2833 

Lacey, C, & Cole, S. 1993, MNRAS, 262, 627. 

Ma, C.-P., & Bertschinger, E. 1995, ApJ, 455, 7 

Mo, H.J., & White, S.D.M. 1996, MNRAS, 282, 347 

Navarro, J., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 



-38- 



Peebles, P.J.E. 1969, ApJ, 155, 393 

Press, W.H., & Schechter, P. 1974, ApJ, 187, 425 

Peebles, P.J.E. 1993, Principles of Physical Cosmology (Princeton: Princeton University 
Press) 

Ricotti, M. 2003, MNRAS, 344, 1237 
Sycr, D., & White, S.D.M. 1998, MNRAS, 293, 337 
Subramanian, K., Cen, R., & Ostriker, J.P. 2000, ApJ, 538, 528 
Xu, G. 1995, ApJS, 98, 355 



This preprint was prepared with the A AS IM^jX macros v5.2. 



